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Motivated by the prospect of Bardeen-Cooper-Schrieffer (BCS) pairing in cold fermionic gases 
we analyze the superfluid phase of 3 fermionic flavors in the attractive Hubbard model. We show 
that there are several low-lying collective pairing modes and investigate their damping due to the 
partially gapless nature of the single-particle spectrum. Furthermore we analyze how these modes 
show up in the density response of the system. Apart from the Anderson-Bogoliubov phase mode 
of the pairing between two flavors, the dynamical structure factor contains signatures of the gapless 
third flavor. This picture is found to be robust against perturbations that break the global SU(3)- 
symmetry of the Hamiltonian. 

I. INTRODUCTION 

The experimental realization of Bosc-Einstcin condensation^ has opened up a new avenue of research in systems 
of ultracold atoms. Recently, improved cooling techniques have also led to degenerate Fermi clouds^'S*^. An active 
search is now under way to implement a BCS transition of degenerate fermionic gases to a superfluid stated, and 
very recently first indications of pairing of ''°K atoms have been reported^iSiS. Although interatomic interactions 
are generally weak, it has been become possible to approach strongly correlated phenomena either via Feshbach 
resonancesSiiSiiiii^ or optical latticesi^. In this work we focus on the latter possibility. In a pathbreaking recent 
experiment, a 3d optical lattice has been used to experimentally realize the transition between a bosonic superfluid 
and a Mott insulator—. In combination with fermionic atoms, optical lattices offer the intriguing perspective of 
studying solid-state phenomena like high-temperature superconductivity in a new contexlji^. 

Compared to electrons in solids, atomic systems offer new internal degrees of freedom that can lead to states of 
matter which do not have obvious counterparts in the physics of interacting electrons. For alkali atoms, nuclear spin 
/ and electron spin S are combined in a hyperfine state with total angular momentum F. While typical electronic 
systems are constrained to SU(2) spin rotational symmetry, the atomic total angular momentum F can be larger than 
1/2, resulting in 2F + 1 hyperfine states differing by their azimuthal quantum number mp. In a magnetic trap, only a 
subset of these 2F+ 1 states (the low-field seekers) can be trapped, but this constraint can be avoided using all-optical 
trapsi^. In fact, coexistence of the three hyperfine states \F = 9/2, mp = —5/2, —7/2, —9/2) of '^^K in an optical trap 
has already been demonstrated, with tunable interactions due to Feshbach resonances between mp = —5/2/ — 9/2 
and mp = —7/2/ — 9/2, respectively^. A situation with strong attractive interaction between all three components 
can be realized e.g. for the spin polarized states with iris = 1/2 in ^Li where the triplet scattering length a — — 2160ao 
is anomalously largei-''. 

Optical lattices are created by a standing light wave leading to a periodic potential for the atomic motion of the 
form V{x) — Vb cos^(fca;i) where k is the wavevector of the laser, i labels the spatial coordinates and the lattice 
depth Vf) is usually measured in units of the atomic recoil energy En = ?i^fc^/2m. In the following we will consider the 
2D case where i = 1,2. It has been showni^ that the Hubbard model with a local density-density interaction provides 
an excellent description of the low-energy physics. Here we are interested in a situation where fermionic atoms with 
A'' different hyperfine states (further on denoted as "flavors" ) m are loaded into the optical lattice. We thus consider 
a Hubbard Hamiltonian 



H 



Here Ui — JT-i,™ is the total number density of atoms on site i which can be written in terms of the fermionic 
creation and annihilation operators, rii^m = cl^^Ci^rn- The interaction (second term in Eq. ^ is invariant under local 
U(A^) rotations of the N flavors with different m. The hopping term of the atoms between nearest neighbors (ij) 
reduces the invariance of iJ to a global V{N) symmetry. Stripping off the overall U(l) phase factor, we arrive at 
the SU(A^) Hubbard model. In the optical lattice the Hubbard parameters are t — (2/^/??) exp (— 2^^) and 

U = EpiaskyJS/'K^^ where ^ = {Vq/ Er)^^^ . as is the s-wave atomic scattering length. 

The fermionic SU(A'^) Hubbard model on the two-dimensional (2D) square lattice was studied for repulsive inter- 
actions J7 > in the large- A^ limiti^ in the early days of high-Tc superconductivity, mainly as a controllable limit 
connected to the then physically relevant case N = 2. Affleck et alji^ already discussed realizations of SU(4) using 
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FIG. 1: a) Hubbard model for three flavors of fermions. The particles can hop with amplitude t between neighboring sites 
of the optical lattice and interact on the same lattice site with an onsite interaction U. b) Three possible types of onsite 
pairings for three flavors of fermions with attractive interaction. In the mean-fleld theory for the SU(3) symmetric case all 
three types of pairings and superposition thereof have the same ground state energy, c) Density of states (DOS) for the three 
fermionic excitation branches of the SU(3)-paired state (schematic, lattice effects neglected). In the gauge with only A12 7^ 0, 
the unpaired flavor 3 is gap less while the two branches involving flavors 1 and 2 see an energy gap A. 



the nuclear spin of ^^Ne. A generalized SU(A^) model could describe orbitally degenerate electronic states in crystals, 
but it is likely that different overlaps between the orbitals pointing in distinct lattice directions will break the SU(A^) 
invariance. In a previous work^S we have shown that exotic states like staggered flux states and partially incompress- 
ible flavor density wave states may occur in the repulsive 2D Hubbard model for finite N > 2. We also pointed out 
that in the attractive case for N > 2, the onsite Cooper pairing breaks more symmetries than just a global U(l) 
symmetry. This gives rise to new collective modes in the paired state. 

In this paper we focus on the attractive regime U < and the question how the spin-1/2 BCS state will be 
generalized for a number of flavors N > 2. We shall assume weak to intermediate local attractions such that a 
treatment within BCS theory is qualitatively valid. We focus mostly on TV = 3 where we analyze the ground state 
manifold using mean-field theory and determine the collective modes by a generalized random phase approximation 
(RPA). We show that the density response obtained by Bragg scattering provides a clear experimental signature for this 
new superfluid ground state. Although our considerations take place in the framework of the Hubbard lattice model, 
the qualitative features of the excitation spectrum should remain the same for continuum systems with attractive 
interactions between three flavors. 



II. MEAN-FIELD THEORY 



Let us briefly describe the results of BCS mean field theory for the SU(A^)-model with attractive contact (Hubbard) 
interaction. Within the mean-field approximation these results are quite general and hold independent of the spatial 
dimension, also for continuum systems without underlying lattice. Parts of the analysis are already contained in a 
brief form in Ref. I20I 

The Hamiltonian we study is 

k,a k.k',q 

a/] 

with a kinetic energy — e_^, flavor indices a, f3 1 . . . N and chemical potentials Ha- V denotes the number of 
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lattice sites. We introduce a pairing mean-field 

k 

and consider even-parity pairing A^p = —Apa (see also Fig. Non-trivial odd parity solutions are not possible for 
local attractions. The mean-field Hamiltonian then becomes 

Hum- = E ^^-fc/j^/^" + 

k,a,0 

For N > 2 the even parity gap functions A^^ transform non- trivially under SU(7V). This means that depending on 
the global gauge of the fermions, Aap takes different values. This is different from the SU(2) case. There even parity 
gap functions arc singlets and invariant under spin rotations and transform with the trivial representation. From a 
particle physics perspective we can compare the BCS pairing in the SU(3) case with the pairing of two quarks or two 
antiquarks. The quarks transform according to the color-SU(3) which is a local gauge symmetry in this case. Bound 
states of two quarks or two antiquarks carry color-SU(3) charge and do apparently not exist in nature. In contrast 
with that pairs of quarks and anitquarks, e.g. pions, can be color-SU(3) singlets. 

Without much group theoretical input we can read off the dimension of the representation in the SU(Af) case: 
an antisymmetric complex N x N matrix Aap has N{N — l)/2 independent complex components. For example for 
iV = 3, Aai3 transforms according to a 3D representation. This is consistent with the decomposition of the product 
3(8)3 = 3®6. Here 3 denotes the irreducible representation under which the flavor spinor [cj: j^, 2' '^k 3) transformed, 
and 3 is the complex conjugate representation of 3. This can be checked explicitely when we write (dropping the 
wavevector indices for a moment) 

1 / 

Aa = -eaPjicpC^) ^ -Ai3 

\ A12 

and analyze the transformation Aq — D{U)af3Ai3 when the fermions are transformed as Ca = Uapcp, where U is an 
SU(3) matrix. Expanding Uap = Sap + iuaf) and Daf3 ~ 5ap + idap for an infinitesimal transformation and using the 
traceless hermitian nature of Uajs and da/s, we obtain dafi = — Wq^, i-e. D = U* . This transformation property implies 
that any 3-vector A„ with lA^p = Aq can by a global gauge change be mapped onto (0, 0, Aq), where A12 — Aq 
and A13 = A23 — 0. In this gauge the mean-field solution is simply the solution for the SU(2) BCS problem in flavors 
1 and 2 plus the third flavor that remains unpaired and thus ungapped. Therefore in the ground state and weak 
coupling iV_F|C/| < 1 

Ao = W^e-i/(^-l^l), 

where Np is the density of states at the Fermi level and W the bandwidth within which the attraction /7 < is active. 
For 3/8 band filling and U = —it, the mean-field Tc is ~ 0.17t. The quasiparticle spectrum contains three branches, 

El -2 = ±^ei + Al and E3 ^ ej: , 

and is independent of the global gauge. Therefore the paired mean-field ground state of the SU(3) BCS problem has 
a full Fermi surface and only two thirds of the excitations are gapped (see also Fig. ^ . 

For SU(2), the ground state is degenerate with respect to the global phase of the gap function, and long- wavelength 
variations of the latter are gapless in absence of long-range forces. In the SU(iV) case we find a higher degeneracy of 
the ground state and more gapless modes as more generators are broken. As argued above, for SU(3) all gap functions 
with the same Aq = J^ap |Aa/3p are degenerate and have the same total density of states. Apart from the global 
phase there are four additional gapless modes, two associated with the internal phases between A12, A13 and A23, 
and two modes modulating IA12I, IA13I and IA23I with fixed Aq. 

The gauge with only A12 7^ and A23 = A13 = makes the symmetry breaking pattern obvious. The original 
symmetry group of the problem SU(3) (g) U(l) has nine generators. This gets broken down to an SU(2) symmetry in 
flavor 1 and 2, leaving A12 invariant, and an additional U(l) that acts on the phase of the unpaired flavor 3. This 
leaves 5 generators broken, yielding the collective modes described above. 

The coexistence of a full Fermi surface with a superconductor should have interesting consequences. Some of these 
will be analyzed in the following sections. 
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Let us briefly comment on perturbations that break the SU(3) symmetry of the Hamiltonian. A weak magnetic 
field will cause a Zeeman splitting of the hyperfine levels with different mp. This is equivalent to having different 
chemical potentials for the three flavors. As the flavor densities become unequal, the mean-field solution for zero 
total momentum pairing retains a nonzero Ao if the Zeeman splitting between two neighboring hyperfine levels m^?, 
mi? + 1 is smaller than Aq. We do not consider the possibility of pairing with non-zero total momentum, leading to 
a so-called Larkin-Ovchinikov-Fulde-Ferrell (LOFF) state^i. For the continuum model it has been foundSS that zero- 
total-momentum pairing gives the lower energy for sufficiently strong attraction. If the interaction remains the same 
between all flavors, the zero-total- momentum pairing will choose those hyperflne levels that have the largest density 
of states Pa near the respective Fermi levels, i.e. if pi > P2 > Ps, the pairing amplitudes pick the solution IA12I = Ag 
and A23 — A13 — 0. In the case with unequal strengths of the attraction between the flavors, the pairing amplitude 
chooses the strongest interaction. This means for —U12 > — f/23 and — C/12 > —U13, we again have IA12I = Aq and 
A23 = Ai3 = 0. 

For SU(4), the antisymetric matrix Aaf3 has 6 independent components A^, i = 1, . . . , 6, which transform under 
a 6-dimensional representation. Then with an SU(4) rotation we cannot always map a general 6-vector A; onto a 
vector (X Aq • (0, 0, 0, 0, 0, 1). Correspondingly the degeneracy of the ground state is subject to more constraints than 
just constant Aq. The mean-field solutions have IA12I = IA34I, jAiaj = IA24I and |Ai4| = IA23I. The single particle 
spectrum is fully gapped. 

III. GENERALIZED RPA FOR THE SUPERFLUID STATE 

Here we describe a generalized random phase approximation (RPA) scheme that allows us to analyze pairing and 
density fluctuations in the attractive SU(3) Hubbard model. 

The theoretical description of the interference of different collective fluctuations of a many-fermion system is a 
complex problem. In principle the interacting two-particle scattering vertex contains all the information about pos- 
sible collective modes. However, for general situations where charge, spin and pairing channels are important, the 
calculation of this object, e.g. using a Bethe-Salpeter equation, is rather involved even if self-energy effects are ne- 
glected, because both particle-particle and particle-hole diagrams have to be taken into account. Fortunately, in many 
situations one can with good confidence restrict the analysis on certain channels and apply a Hubbard-Stratonovitch 
(HS) decoupling of the interaction in physically reasonable way. In order to motivate our approach let us sketch our 
expectations regarding the collective excitations in the SU(3) case with local attractive interactions. 

The pairing channel contains soft collective modes due to the spontaneous symmetry breaking and we expect that 
the density of the fermions will couple to these modes. For example for superconductors with long range Coulomb 
interactions it is knownS^ that the density response "eats" the Goldstone phase mode of the superconductor by 
pushing it up to the plasma frequency, which is the collective mode of the charge density. In our case of an attractive 
Hubbard interaction there is no long range force, but the coupling of the phase mode to the partially gapless density 
fluctuations could still lead to observable effects. In particular one may wonder whether the phase mode gets strongly 
damped by the third flavor. A similar obliteration of a collective mode occurs in the normal state for U < 0. There 
in RPA the zero sound mode of the three flavors is located inside the particle-hole continuum and thus it is strongly 
Landau damped. In the case of three flavors with pairing of two flavors, the ungapped flavor could still maintain a 
zero sound mode that sharpens below the transition as a part the quasiparticles become gapped and the damping is 
reduced. 

In the flavor channel, in analogy with the spin channel of the SU(2) case, the sign of the effective interaction is 
reversed and for U < there are two flavor modes with linear dispersion above the particle-hole continuum. In the 
normal state these modes can be understood as relative density oscillations oc pi + p2 ~ 2p3 and cx pi — p2 of the three 
flavors that keep the total density pi + p2 + Ps fixed. Hence these flavor modes are invisible in the density response. 
In the superfluid state with only A12 7^ however, the density response of the gapped flavors 1 and 2 is different from 
that of the gapless flavor 3, and the mode (x pi + p2 — 2/33 couples to the total density. This could make the flavor 
mode observable in the density response. 

Summarizing this short discussion, we state that a Hubbard-Stratonovitch decoupling scheme should include the 
pairing degrees of freedom and the individual densities of the three flavors. Let us now describe a formalism that 
allows a treatment of the issues raised above. It is a straightforward generalization of the methods described e.g. in 
Ref. 0. We start with the attractive Hubbard Hamiltonian of Eq. |21 Next we promote the fermionic operators to 
Grassmann numbers and introduce the combined index k = {iuj,k). Then the partition function for the system is 
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given by 



/ 



DcDc exp 



U 



^ Cfca(ia; + eg_^ - /la) Cfea - — ^ Ck+q,aC-k,0C-k',/3Ck'+q 



k.a 



(3) 



Now we decouple the local attraction in pairing, density and flavor channels. We write 



(4) 



Here the first term is decoupled in the onsite pairing channel with a coupling constant Up <0, while the second term 
with Up ^ U — Up < is decoupled with local Hubbard-Stratonovitch fields coupling to the densities of the three 
fiavors. To this end we write in the pairing channel 



exp 



n(}. 



2\U\ 



/3,i| 



/ 



DAiDA* exp [-Sa{c, c, A, A*)] 



(5) 



(6) 



Searching the saddle point with respect to the Hubbard-Stratonovitch field Aaj3 reproduces mean-field theory. Let us 
write 



Aa/3,i = Aai3 + Aai3^i 



(7) 



where Aap denotes the static site-independent saddle point solution and Aafi.i are the fluctuations around it. The 
saddle point solution Aa/s can be absorbed into the mean-fleld action Smf- Next we split up Aajs in real and imaginary 
parts, Ac/3,i = AJj^ j + zA^^ j. Going over to wavevector-frequency space we obtain the pair fluctuation action 

5a (c, c, A^ A') = E h+<i,aC-kAK0{'l) + iKpil)) + (A^X?) - «Aj:^(g))c_fc,^Cfe+,,„] 



+ ^E[|A;M9)l' + |Aa/3(9)l'] • 



(8) 



In the flavor and density channel we write the interaction as a sum of squares, 



V3 



ni.i + n2.i - 2n3,, 



12 



■111,,. - "2.1 



V4 



and use the HS identity 
exp I \Up 



niA + n2,i + nz,; 



niA - n2,i 



= j Dpt,i Dpi23,i Dpi2,i |-^|^ [Pt,i + P?23,i + Pl2,i] 



niA + n2A + ns i . nn + n2,i - 2n3 , . ni.i ^ ^2,; 

+ Pt,i 7= \-ipl23,i 7= \-iPl2,i 



I 



V3 ■ V12 

Dpt,i Dpi23,i Dpi2,i exp [-Sp{c, c, pt, p\2, P123)] 



V4 



(9) 



with the real site-dependent HS fields pt^i, pi23,i and pi2,i- Here pt,i couples to the total density, and pi2,i and P123 
couple to out-of-phase oscillations of the flavors. We will denote these as flavor modes. For all situations we study, 
the static components of the flavor fields and also the density field pt{q, a; = 0) away from q^O remain zero, i.e. we 
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do not study flavor or charge density wave states. In the case where Up — > 0, the internal pa fluctuations get pinned 
at zero and do not contribute to the partition function. In Fourier space we have 



Sp{c, c, pt,pi2,pi23) = - ^Pa{q)Na{-q) 



1 



4|C/p| ^ 



\Pa{q)\' 



(10) 



where the index a labels the channels t, 12 and 123 and Nt{q) = vJ2k a'^k+q,aCk,a, -^12(9) = 
W T^k,a (cfe+g,iCfe,i - Ck+q,2Ck,2) and iVi23(g) = -7=r7 J2k,a (ck+q.iCk,! + Ck+q.2Ck.2 " 2ck+q,3Ck,3) ■ In the functional 



integrals over the HS flelds A^J^^ and pa we have to keep in mind that their real space versions are real fields. Therefore 

we have Al/f^'*iq) = Kipi-l) and pl{q) = Pai-q)- 

Finally we add an external total density field p{q). It couples to the total density of the fermions in the same way 
as pt- However as external field it is not integrated over in the partition function and there is no contribution oc 1/Up. 
The action for this coupling is 



S'cxt. (c,c) = -^^p{~q) 



(11) 



The quadratic fermionic parts of these four contributions to the action, 5'mf, Sa, Sp and Sc^t. can be cast into a 6 x 6 
Nambu matrix form using a 6- vector = {ck^i, Ck^2: c^^a, c_fci, c_fc2: C-^a). This way we have doubled the fermionic 
degrees of freedom and have to restrict the /c-space sums to half of the Brillouin zone. Up to the non-fermionic parts 
the action can be written as (using the summation convention in wavevectors fc, the Nambu-flavor index a and the 
fluctuation index A) 



QC f^k/y^—l\kk /^k 



Ck jpkk' A rkk' /^k' 



(12) 



Here G ^ is the inverse Green's function of the saddle point theory, containing A^p. "^f^p is the vertex factor 
coupling a flavor or pairing HS field Fa with wavevector /frequency fc — fc' = q to the fermion bilinear. As we have 6 

pairing fields A^^*, 3 flavor/density HS fields and 1 external source field, A = 1, . . . , 10. Integrating out the fermions 
yields the determinant of this quadratic kernel, or its logarithm if we put the result back into the exponent. Then we 
use log det Sp = tr log Sp and obtain the following contribution to the action for the fluctuations 



5|, = -trlog {G-'ti;-Ff 



(13) 



Our goal is to expand this expression to second order in A^^ . Together with the contributions from the last term in 
the action (jSJ this will give a quadratic action for the fluctuations. The terms linear in the fluctuations cancel provided 

we expand around the saddle point. Next we write G^^ —Fa f ~ G^^{1 — GFa ^ f) and use log(l — x) — — x— ^ — . . .. 
Thus in order to extract the fermionic contribution to the fluctuation action, we have to evaluate 



S 



r'k Tpq A rq r^k-q Tp~q A' 
'^al3"^A J I3"P''^!3'P"'^ A' J f3- 



J 3"'c 



(14) 



Here we have used the fact that G^'^, is diagonal in k, k' and denoted k—k' = q. Diagrammatically, these contributions 
correspond to one-loop bubbles made out of two Green's functions and perturbations F^ at the vertices. Including 
the non-fermionic contributions the quadratic fluctuation action reads 



S 



(2) 



n 



d-A^AA' + ^fl'"aG^p"G^,3t^. 



A' f-q 
Jf3"l3' 



Fa^ 



(15) 



Here the diagonal factors dA are (La — V"/|C/p| for the pairing HS fields and cLa = ^/4|?7p| for the flavor-density HS 
fields. For the external source field pifi)-, dA=io = 0, and the only contributions are due to the fermionic bubbles. It 
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is convenient to group the HS fields into a 10- vector 



Hi) 



With this Eq. [Threads 



Al2(g) 
Ala (9) 
Ai3(<?) 

pt{q) 
Pi23iq) 

p{q) 



(16) 



(17) 



This quadratic action couples the internal degrees of freedom among each other and to the external field. Like in any 
other RPA, only modes with the same frequency /wavevector q interact with each other. If we are interested in the 
Green's function of a specific fluctuation F like the amplitude or phase of a Aq,^ in absence of the external field, we 
integrate out the remaining 8 internal degrees of freedom yielding 



M 



FF 



M 



FF 



E 



Mff'M 



F'F' 



F"F ■ 



(18) 



F" = l, 



Here the sums F' and F" run over the other 8 internal degrees of freedom, excluding F and the external source p. 
The Green's function for fluctuation F is given by the inverse of Mpp{q), and the corresponding spectral function is 
obtained after analytical continuation ^ uj + id as Ap{q,uj) = Im [Mpp{q,u!)] ^. 

Regarding experimental manifestation of the collective excitations we are interested in the total density response. 
This is given by the renormalized quadratic part in p{q) after integrating out the 9 internal HS fluctuation fields. 
Doing the Gaussian integral yields 



E 

F = l,..-9 
F' = l,...9 



MpF'M 



F'F' 



.Mf 



"P ' 



(19) 



where this time the indices F' and F" go through the 9 internal degrees of freedom. The total density susceptibility 
Xp{q) is then given by the double derivative of the remaining exponential with respect to the source fields p{q) and 
p{—q) and the sources set to zero afterwards. This simply yields 

xp{q) = M;^(g) . 

The imaginary part is the total density spectral function Im Xpiq^^)- 

The density susceptibility Xp(g, w) is relevant for experiments in ultracold atoms, where the dynamical structure 
factor 



5(g,a;) = El(/'l)„ol''5(^-^"o) 



(20) 



can be measured via Bragg scattering^^'^. Here ojno = En — Eq denotes the many-body excitation energies relative 
to the ground state. Within this technique, two laser beams with wavevector- and frequency-difference k and uj are 
used to create a time-dependent potential {light grating) 

Vn,od{x,t)^Vcosiqr-ujt) (21) 

leading to scattering of atoms due to two-photon processes. The resulting transition rate per atom is given by 



W/N = -V^S{q,Lu). 



(22) 
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At T = the structure factor is related to the density susceptibihty according to 

ImXp(g,c^) = ~TT{S{q,uj) ~ (23) 

i. e. at positive frequencies we have S{q,uj) — —^lmxp{q,uj). From our calculation of the density response in Eq. H19|) 
we can thus extract the Bragg scattering rate. 

In the decomposition Q we will later choose Up ~ Up = U/2, but in principle this decomposition is arbitrary and 
all choices are equally exact as long as the full action is treated. Yet within the Gaussian approximation we have to 
resort to in the further evaluation of the partition function, different decompositions can in principle lead to different 
saddle point solutions. This is a well known ambiguity of the Hubbard-Stratonovitch formalism. In our case however 
the qualitative results do not depend on the precise partitioning of the interaction. Sufficiently away from half band 
filling, the ground state in the saddle point approximation will always be given by the BCS mean-field solution, 
provided Up < 0. The gap value Aq will be reduced corresponding to a weaker pairing interaction Up = U — Up. 

IV. COLLECTIVE MODES FROM PAIRING FLUCTUATIONS 

To begin with, let us set Up = and decouple the interaction in the pairing channel only. This allows us to 
analyze the pairing fluctuation as internal degrees of freedom without interaction to other collective fluctuations. The 
coupling to the flavor densities and the total density response will then be analyzed in the following section. Unless 
stated otherwise, all numerical results shown in the following section were obtained for the parameters Up — — 4i and 
(n) ~ 3 • 0.185 — 0.556 particles per site. The mean-field critical temperature for this case is Tc « 0.07t and Aq ~ O.lt 
at zero temperature. We analyze the spectral functions of the fluctuations for wavevectors q — {q, 0) with small q 
along the x-axis. As the Fermi surface is anisotropic, different directions in (f-space give slightly different results. 
However we have checked that the qualitative features we show in the following are common to all cuts through the 
Brillouin zone. 



A. SU(3)-symmetric case 

First we consider the SU(3) symmetric case. As explained in Sec. II, the mean-field ground state is degenerate for 
all choices of saddle points Aap with -i p l^a/^P = ^o- For simplicity we will always choose the gauge with only 
Ai2 = — A21 = Aq with real Aq and A23 = A13 = 0. Then the A-sector of the matrix M in Eq. El is diagonal and 
we can just read off the spectral functions for the amplitude mode /S.\2 the phase mode A22 by inverting the 
corresponding matrix entries and taking their imaginary parts. For A13 and A23 the saddle points are at the origin. 
Therefore their real and imaginary parts fluctuate identically and have the same spectral function. 

Due to the constraint i ^ l^a/3p = ^0 '^^ expect soft modes for the 5 possible variations of A^^g that leave Ag 
unchanged. Since only A12 7^ 0, these modes are the phase of A12, and real and imaginary fluctuations of A13 and 
A23. In addition to that there will be an amplitude mode of A12 at twice the gap frequency. 

The results of the generalized RPA are shown in Fig. |21 In the left panel we find the phase mode of flavors 1 and 2 
and the gapped amplitude mode. The Goldstone phase mode is gapless and start with a linear dispersion at small q. 
The velocity c of the mode for our lattice model is very close to the continuum result c — vp/^/^ in two dimensions, 
where vp is the Fermi velocity along q {c = vp/V^ in 3D)^'^. The amplitude mode has a gap of 2Ao and is not sharp 
as it can decay into quasiparticle pairs above the gap edge. In this approximation the phase mode is undamped for 
uj < 2Aq. However we will show below that remaining density-density interactions for Up ^ Q lead to a broadening of 
the phase mode. We also observe that the soft modes of A 13 and A23 get strongly broadened above w = Aq. Flavor 3 
is ungapped and therefore these modes can excite quasiparticle pairs of flavor 1 and 3 (or 2 and 3) already above Aq. 

B. Broken SU(3)-symmetry 

Here we briefly describe how two types of SU(3)-symmetry breaking in the Hubbard Hamiltonian affect the internal 
pairing modes. 

The first case we consider is unequal values of the attraction between the flavors, for example — U12 > —U13 > — 1/23. 
Then the mean-field solution chooses to pair the two flavors with the strongest mutual attraction. In the speciflc 
example this corresponds to A12 7^ 0, while the other two pairing amplitudes are zero, A13 — A23 = 0. 

An analysis of the internal pairing modes gives the following picture. The modes of this dominant Ai2-pairing 
channel remain unchanged. In particular the phase mode remains gapless. This is plausible as the unequal strength of 
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FIG. 2: Collective pairing modes in the SU(3)-symmetric case with onsite s-wave pairing A12 = Ao, A13 = A23 = for 
Up = — 4t, filling n — 0.556 and T — O.Olt. The steps in the data are due to the numerical evaluation for a finite number 
of wavevectors and frequencies. Left plot: Spectral weight corresponding to fluctuations of A12. The soft linear mode is the 
Anderson-Bogliubov phase mode, the weaker gapped feature at 2Ao is due to amplitude fluctuations. Right plot: Spectral 
weight of A13 or A23 fluctuations. Since the mean fields A13 = A23 = 0, real and imaginary parts have the same fluctuation 
spectrum. Landau-damping due to the gap less fiavor 3 starts at cu — Aq. 
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FIG. 3: Collective pairing modes in the SU(3)-broken case with onsite s-wave pairing and A12 = Ao, A13 — A23 = 0, U12 ~ 4t, 
U13 = 3.6t, and U23 = 2t, T — O.Olt. Left plot: Spectral weight corresponding to fiuctuations of A13. The dark feature is due 
to two gapped modes corresponding to fiuctuations of the real and imaginary part of A13 just below the particle-hole (PH) 
continuum. Right plot: Spectral weight of A23 fluctuations. The attraction U23 = 2t in this channel is not strong enough to 
pull the modes below the PH continuum and the modes are almost wiped out by the PH-pairs. 



the interactions does not lift the degeneracy of the ground state with respect to the phase of A12. If the attraction in 
a subdominant channel is strong enough, the collective oscillations of the pairing field in this channel get pulled below 
the particle-hole continuum and hence get sharp. This is shown in Fig. 13 where the attraction in the 13-channel is 
strong enough to yield a sharp mode below Aq, while the mode in the 23-channel with the least attractive interaction 
remains above the Ag-threshold and is strongly damped. It may actually be interesting to see if the fluctuation energy 
(i.e. the frequency) of these modes could be lowered by a modification of the mean-field state. Switching on a nonzero 
pairing amplitude in the subdominant channel will have an opposite effect as it pushes the modes further up. 

Next we consider the effects of a Zeeman splitting between the hyperfine levels. This corresponds to different 
chemical potentials or, equivalently, different densities of the three fiavors. For simplicity let us assume a Zeeman 
shift +Ez for flavor 1 and —Ez for flavor 3, respectively, while the energy of fiavor 2 remains unchanged. For band 
filling of less than half filling, the pairing will select A12 = Aq, as then flavors 1 and 2 have the larger densities of 
states. As mentioned in Sec. II, the mean-field solution breaks down only when the band splitting Ez reaches the 
energy gap of the pairing Aq. Even in this case of SU(3) symmetry breaking, the degeneracy of the mean-field state 
with respect to the phase of A12 is not lifted. Therefore the phase mode of A12 remains a Goldstone mode. As the 
symmetry between the flavors is broken, the modes in the unpaired channels acquire gaps of the order of the Zeeman 
splitting. 
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FIG. 4: Density response for the superfluid state at T = O.Oli and A12 — Aq, A13 = A23 = in the case without remnant 
density-density interactions, Up = 0. Left plot: Im x''('7i '^)- The upper hnear feature is the upper boundary of the particle- hole 
continuum of flavor 3, the lower feature is the density signature of the Anderson-Bogliubov phase mode. The amplitude mode 
and the modes involving A23 and A13 do not appear in the density response. Right plot: Cut through Im x''('7j for q/o, = 0.08 
(see mark in left plot). The solid line is the density response of flavor 1 and 2, exhibiting the phase mode. The dashed line is 
the non-interacting particle-hole continuum of flavor 3. 

V. DENSITY RESPONSE 

Having studied the internal pairing modes of the SU(3) superfluid, we now ask how these excitations couple to 
density oscillations. This is important in order to find experimental signatures of the collective excitations via Bragg 
scattering, as explained in section UTTI In fact, for the SU(2) case the appearance of the Anderson-Bogoliubov phase 
mode has been suggested as a promising way to detect the onset of fermionic superfluidity^^ . 

A. Pure collective pairing fluctuations in the total density response 

First, in order to separate the different effects, we study the density response when the interaction is decoupled 
in the pairing channel only, i.e. we keep Up = Q in Eq. 01 This means that the only collective excitations stem 
from pairing fluctuations. Analyzing the imaginary part of Xp((f, a;) we find that only one of the six collective pairing 
modes show up in the density response. This mode is the phase mode of the nonzero pairing amplitude. This is 
not surprising. If the pairing amplitude in a certain channel is zero (say A23 = 0), there is no superfluid density 
that could react to a phase fluctuation. Thus there is no quadratic coupling between density and phase. The second 

observation is that the density signature of the phase mode in the density response goes to zero for g — > 0. This can 

2 

be understood as follows. From the continuity equation we get 5p oc ^S(j), i.e. if a; ^ 5 the density change Sp induced 
by a phase variation 6(j) vanishes like q. Consistent with that we find that the density response of the phase mode 
becomes smaller for g ^ 0. The amplitude mode of the paired flavors a,t uj — 2Ao is not reflected in the density 
response, just as in the case of only two flavors^!. 

Thus for the case with only A12 ^ the density response consists of two features. The first is a sharp phase mode 
that - in RPA - is undamped below 2Ao, and the second is a particle-hole continuum of the 3rd unpaired flavor. 
Notice that in 2D the upper edge of the continuum is rather sharp. This continuum feature distinguishes the density 
response of the 3-flavor case from that of the conventional 2-flavor superfluid. Below we will see how remaining 
density-density interactions affect this feature. 

B. Additional density-density interactions 

Next we want to analyze the effects of remaining density-density interactions. So far we have decoupled the 
interaction in the pairing channel only, but as argued above collective density fluctuations can play a role as well. For 
that purpose we set Up = Up = —At. 

In the normal state, we find two degenerate flavor modes corresponding to P123 and pi2 which lie somewhat above 
the particle-hole continuum and are therefore undamped at low energies. However these modes do not show up in 
the density response, pt does not exhibit a sharp mode, as can be expected for attractive interactions. In this case. 
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FIG. 5: Density response for the normal state at T = 0.2t ~ 3Tc and A12 — A13 = A23 — 0, for remnant density- density 
interactions Up — Up — —At. Left plot: Im x''{q,uj) in the normal state when density-density interactions are included. Right 
plot: Cut through Im x''[q, ijj) for q/a = 0.03 (solid lines) and q/a = 0.08 (dashed lines, see marks in left plot). The thicker lines 
correspond to the interacting case and the thin lines with the sharper peak show the non-interacting particle-hole continuum. 




FIG. 6: Density response for the superfluid state at T = Q.Qlt and A12 = Aq, A13 = A23 = 0, again for Up — —At. Left plot: 
Im x'^i'ii^)- The upper linear feature is the signature of the flavor mode of P123, the lower feature is the density signature 
of the Anderson-Bogliubov phase mode. The amplitude mode and the modes involving A23 and A13 do not show up in the 
density response. Right plot: Cut through Im x''{q,uj) for q/a — 0.03 (solid line) and q/a — 0.08 (see marks in left plot). The 
inset is a zoom-up of the shaded region around ui/t = 0.16. The dashed line is the density response at q/a = 0.08 and the 
crosses coinciding with it are the data points for the rescaled spectral function of the flavor mode P123. For comparison, the 
other solid curve with the peak pointing downwards is proportional to the non-interacting particle-hole continuum. 



instead of creating a zero sound mode, the density-density interactions smear out the sharp upper edge of the bare 
particle-hole continuum, and the response of the total density exhibits only a broad peak. 

In the paired state, the phase mode of A12 is the dominating feature in the total density response. In comparison 
to the results without density-density interactions we find that the coupling to the density fluctuations causes slight 
damping of the mode. This can be seen by comparing the spectral functions in Figs. ^ and El For our parameters 
and Up — Up the phase mode is still clearly visible and distinguishes the density response of the paired state from the 
normal state result. In principle the temperature dependence of the width of this mode provides a way to detect the 
gapless third flavor. In addition to that, the density-density interactions renormalize the mode frequency to somewhat 
lower energies. In other words, in the paired state the gapping of flavors 1 and 2 sharpens up the zero sound mode that 
is pulled down by the interactions into the particle-hole continuum of flavor 3. Apart from the phase mode, there is a 
second feature in the density response that resembles the upper edge of the non-interacting particle-hole continuum. 
A closer comparison shows however that this second peak lies a bit above the free particle-hole continuum and is due 
to the flavor mode in the spectrum of pi23- This mode involves density variations (x ni +n2 — 2713. Since for A12 7^ 
flavor 3 has a distinct spectrum from flavors 1 and 2, this mode now couples to the total density ni -I- 77.2 -I- 713. Without 
the P123 mode the density response would apart from the phase mode only exhibit a broad peak at the upper edge of 
the particle-hole continuum analogous to the normal state. 

The data described above were obtained for equal partitioning of the interaction with Up = Up. We have checked 
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that other sphttings of the local attraction with Up ^ Up lead to the same qualitative results. For example reducing 
Up simply interpolates between the structure factor shown in Fig. 0]and the data shown in Fig. |H| The damping of 
the phase mode increases with \Up\. The additional flavor mode of P123 merges continuously into the particle- hole 
continuum of flavor 3 when we let Up go to zero. Larger values of \Up\ pushes the mode further above the continuum 
by increasing its linear dispersion with q. Therefore the damping of the Anderson-Bogoliubov mode as well as the 
/5i23-feature near the upper edge of the particle-hole continuum of flavor 3 are robust results that do not depend on 
the ambiguity of the HS transformation in several channels. 

We close the section by commenting on perturbations that break the SU(3)-symmetry of the Hamiltonian, like 
unequal interaction strengths between the flavors and Zeeman splitting. Their effect on the pairing modes was 
considered in the previous section. As the phase mode remains robust with respect to these perturbations and 
the single-particle excitations do not become altered qualitatively, we do not find significant changes in the density 
response as long as the mean-field state remains unchanged. 



VI. CONCLUSIONS 



We have analyzed pairing and density fluctuations in the BCS paired state of the attractive SU(3) Hubbard model. 
This model can potentially be realized by trapping ultracold fermionic atoms with 3 different hyperfine states (flavors) 
in an optical lattice. We have studied the specific case of the two-dimensional square lattice, but most results also 
hold for 3 spatial dimensions and even in the absence of an optical lattice (pure harmonic trapping). 

Even parity onsite BCS pairing in the SU(3) Hubbard model leads to a partially gapped quasiparticle spectrum. 
In the gauge where only A12 ^ 0, two flavors get gapped and the third flavor retains its full Fermi surface. As the 
pairing amplitude transforms non- trivially under 811(3)® U(l), in addition to the usual Anderson-Bogoliubov 
phase mode there are 4 other Goldstone modes. However only the phase mode corresponding to the two gapped 
quasiparticle branches couples to the density of the fermions. Hence only one of the 5 Goldstone modes is visible in 
the dynamic structure factor 5((f, w). This should be observable e.g. in Bragg scattering experiments. For the case of 
two flavors with local attractions, it has been shown that the emergence of the phase mode provides a clear signature 
of the pairing. This occurs because in the normal state the zero sound mode gets wiped out by low-lying particle-hole 
excitations which are gapped in the paired state. If a third flavor is present this mode remains broadened at all 
temperatures due to the damping by the ungapped quasiparticle branch. However, it is clearly visible in the Bragg 
response. 

In addition to the phase mode we find another marked feature in the density response which arises due to a flavor 
mode where the densities of the two gapped flavors oscillate out of phase with the density of the ungapped flavor. 
This mode has a linear dispersion and lies somewhat above the particle-hole continuum of the gapless branch of the 
fermionic excitation spectrum. The separation from the continuum depends on the strength of the remnant density- 
density interaction in the paired state. Together with the visible damping of the phase mode at sub-gap frequencies 
at all temperatures, this flavor mode feature should render a clear signature of the ungapped part of the quasiparticle 
spectrum and thus of the presence of three degenerate flavors in the system. 

We have furthermore investigated the effect of perturbations that break the SU(3) symmetry of the Hamiltonian. 
For example a Zeeman splitting, leading to different chemical potentials of the three flavors, only destroys the pairing 
when it becomes comparable to the energy gap due to pairing, Aq. The density response remains qualitatively the 
same as long as the mean-field state is unchanged. This gives confidence that the described phenomena should remain 
robust over a wider parameter range. 

Acknowledgments: We thank W. Ketterle, P.A. Lee, M. Salmhofer, R. Zeyher, P. ZoUer and M. Zwierlein for useful 
discussions. W.H. was supported by a Pappalardo fellowship. 



^ For a review see special issue of Nature, 416, 206 (2002). 

^ B. DeMarco and D.S. Jin, Science 285, 1703 (1999). 

^ A.G. Truscott et ai, Science 291, 2570 (2001). 

* F. Schreck et ai, Pliys. Rev. Lett. 87, 80403 (2001). 

^ Z. Hadzibabic et al. Phys. Rev. Lett. 88, 160401 (2002). 

® H.T.C. Stoof et ai, Phys. Rev. Lett. 76, 10 (1996). 

C. A. Regal, M. Greiner, and D. S. Jin, Phys. Rev. Lett. 92, 040403 (2004). 

* M.W. Zwierlein et al., Phys. Rev. Lett. 91, 250401 (2003). 

^ C.A. Regal and D.S. Jin, Phys. Rev. Lett. 90, 230404 (2003). 

" S. Inouye et ai, Nature 392, 151 (1998). 



13 



^ E. Timmermans et al, Phys. Lett. A 285, 228 (2001). 
^ M. Holland et al, Phys. Rev. Lett. 87, 120406 (2001). 
^ D. Jaksch et al, Phys. Rev. Lett. 81, 3108 (1998). 
" M. Greiner et al. Nature (London) 415, 39 (2002). 
^ W. Hofstetter et al, Phys. Rev. Lett. 89, 220407 (2002). 
^ S.R. Granade et al, Phys. Rev. Lett. 88, 120405 (2002). 

E.R.I. Abraham et al, Phys. Rev. A 55, R3299 (1997). 
* I. Affleck and J.B. Marston, Phys. Rev. B 37, 11538 (1988); J.B. Marston and I. Affleck, Phys. Rev. B 39, 11538 (1989). 
^ I. Affleck et al, Nuc. Phys. B366, 467 (1991). 

° C. Honerkamp and W. Hofstetter, cond-mat/0309374 (to be published in PRL). 
A. I. Larkin and Yu. N. Ovchinnikov, Sov. Phys. JETP 20, 762 (1965); P. Fulde and R. A. Ferrell, Phys. Rev. 135, A550 
(1964). 

'2 W.V. Liu and F. Wilczek, Phys. Rev. Lett 90, 047002 (2003). 

N.N. Bogoliubov, V.V. Tolmachov, and D. V. Sirkov, A New Method in the Theory of Superconductivity, Fortschritte der 

Physik, Vol. 6, 851 (1958); P.W. Anderson, Phys. Rev. 112, 1900 (1958). 
^ N. Nagaosa, Quantum Field Theory in Condensed Matter Physics, Texts and Monographs in Physics, Springer- Verlag Berlin 

Heidelberg, 1999. 

J. Stenger et al, Phys. Rev. Lett. 82, 4569 (1999). 
^ D. Stamper-Kurn and W. Ketterle, Les Houches Lecture Notes, Session LXXII (1999). 
P. B. Littlewood and C. M. Varma, Phys. Rev. B26, 4883 (1982). 



